Orbital angular momentum spectrum of Bessel-Gaussian modes, 
as generated by spontaneous parametric down-conversion 
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The Bessel-Gaussian modal spectrum, generated in spontaneous parametric down-conversion of 
a Gaussian pump beam, is considered. This is done by first deriving a general expression for the 
true probability of detecting specific transverse spatial modes for the pairs of photons generated in a 
spontaneous parametric down-conversion process. These expressions are applied for Bessel-Gaussian 
modes in type I phase matching with collinear, degenerate down-converted beams. The result shows 
that a broad orbital angular momentum spectrum can be obtained for Bessel-Gaussian modes under 
certain conditions. 



I. INTRODUCTION 

Although the connection between the spatial modal 
profile of an optical beam and orbital angular momen- 
tum (OAM) was initially made with specific reference 
to Laguerre- Gaussian beams pj, the same property ap- 
plies to any optical beam with a rotationally symmetric 
intensity profile. As a result, Bessel beams H, Q (or 
Bessel-Gaussian beams also have quantized amounts 
of OAM associated with them. Each photon in such a 
beam carries an amount of OAM equal to £H, where I is 
the azimuthal index of the mode. 

The use of the OAM eigenstates of photons in quan- 
tum information science became attractive after it was 
shown, theoretically 0, @ and experimentally @, that 
OAM is conserved during spontaneous parametric down- 
conversion (SPDC). A consequence of the conservation 
of momentum in the SPDC process, the conservation 
of OAM implies that the azimuthal indices of a pair of 
down-converted photons add up to that of the pump 
beam. As a result the pair of down-converted photons 
are naturally entangled in terms of their azimuthal in- 
dices or OAM. (This is not surprising, because entan- 
glement in any spatial modal basis implies entanglement 
in every other spatial modal basis that are related to 
the former via unitary transformation. Since momentum 
conservation imply entanglement in terms of the plane 
wave basis, the entanglement in terms of OAM eigen- 
states is inevitable.) Presently, SPDC is the preferred 
method to prepare entangled photonic states for quan- 
tum information processing. 

Entanglement is a desirable property in quantum infor- 
mation application. It is used in quantum ghost imaging 
Q , quantum cryptography 0, [l(| and in quantum com- 
puting algorithms [1J|. Pairs of photons can also be en- 
tangled in terms of polarization, using SPDC with type 
II phase matching. In this case the Hilbert space is two- 
dimensional and the photonic quantum states (qubits) 
are all represented on a Bloch sphere. On the other hand, 
any transverse spatial modal basis defines an infinite di- 



mensional Hilbert space, similar to the plane wave basis, 
which it replaces. A large Hilbert space can represent 
more information in quantum information applications. 

Due to the increased information capacity and the fact 
that SPDC naturally produces OAM entanglement, the 
use of OAM eigenstates for quantum information pro- 
cessing became an attractive proposition. Currently the 
Laguerre-Gaussian modes are the popular choice for the- 
oretical analyses of quantum information systems based 
on OAM. Ironically, the modal basis that are actually 
used in quantum information experiments are very sel- 
dom true Laguerre-Gaussian modes. Although, the he- 
lical phase of these modes can be manipulated and de- 
tected using linear optical systems, the manipulation of 
the full Laguerre-Gaussian modes is not without its chal- 
lenges. The processing and detection of these modes are 
made complicated by the fact that they carry another in- 
dex — the radial index — that governs the radial shape 
of the intensity profile. Each azimuthal index is therefore 
not associated with a unique quantum state, but with an 
infinite subspace of the total Hilbert space. If the radial 
part of its mode changes, a photon propagating through 
an optical system may suffer a loss of the quantum en- 
tanglement, even if the azimuthal index is unaffected. 
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FIG. 1: Experimental setup to detect OAM eigenstate after 
SPDC, showing one beam with a helical phase of I = 3. 

A widely used experimental technique for the detection 
of OAM eigenstates is to use a spatial light modulator 
(SLM) that contains the helical phase associated with a 
specific azimuthal index [HI (See Fig. [J). With this SLM 
the beam is (de)modulated with the complex conjugate 
of its helical phase function, thus removing the helical 
phase from the beam. The resulting beam is then coupled 
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into a single mode fibre to extract the pure Gaussian 
modal component. The light inside the fibre is guided 
onto a photon detector to register the probability for the 
observed quantum state in terms of a photon count. 

If the mode prior to the SLM is a pure Laguerre- 
Gaussian mode with a nonzero azimuthal index the over- 
lap with the mode of the single mode fibre, which re- 
sembles a pure Gaussian, is very small, leading to a very 
small coupling efficiency and very small photon count. 
The reason for the low coupling efficiency is that the ra- 
dial profile of a Laguerre-Gaussian mode with a nonzero 
azimuthal index does not resemble a pure Gaussian. For 
this reason most quantum information experiments do 
not employ actual Laguerre-Gaussian modes. Since what 
is coupled into the single mode fibre is very close to a 
pure Gaussian, the mode just behind the SLM that con- 
tributes to the observed photon count must also resemble 
a pure Gaussian. 

Effectively the physical modes that are being used can 
be viewed as a projected subspace of the full Hilbert 
space of Laguerre-Gaussian modes. The azimuthal in- 
dex is retained, but the radial dependence is projected 
onto something that resembles a pure Gaussian. This 
projection operation is the combined effect of the opti- 
cal system and the detection process that is used in the 
quantum experiment. Imperfections of the optical sys- 
tem, such as its finite transfer function, together with 
the effect of the phase singularity of the helical phase on 
the intensity profile, produce uncontrolled and often un- 
predictable effects on the radial part of the mode. The 
radial part of the modal content of the quantum states 
is therefore often poorly defined and the effect thereof 
on the quantum state during quantum processing is of- 
ten simply ignored. The nature of the radial part of 
the mode does not only affect the coupling efficiency and 
therefore the photon count, but also, more importantly, 
the fidelity of the quantum state. 




FIG. 2: Functions for the generation of Bessel- Gaussian 
modes, (a) Slit-ring aperture with a helical phase of I — 3. 
(b) Conical phase function with a helical phase of £ — 3. 

As the complexity of quantum information systems in- 
crease, the requirements for the fidelity of the quantum 
states can be expected to become more demanding. The 
radial dependence of the spatial modes will have to be 
treated with more care. To this end, Bessel-Guassian 



(BG) modes [H provide better control of the radial de- 
pendence. The radial index associated with Laguerre- 
Gaussian modes is replaced by a continuous scaling pa- 
rameter for the radial part of the BG modes. In practice 
the radial part of a BG mode is generated in the Fourier 
domain by an annulus or a ring-slit aperture, as shown 
in Fig. , the radius of which determines the scaling 
parameter. Alternatively, BG modes can be generated 
by a conical phase function, as shown in Fig.[2fb), where 
the implied conical angle determines the scaling param- 
eter. Both these methods are combined with a helical 
phase function to determine the azimuthal index of the 
BG mode. 

The conical phase method [Fig. Hfb)] is more efficient 
because in this case both the azimuthal and radial de- 
pendences of the mode are determined by phase func- 
tions, which do not introduce a loss in optical power. 
In contrast, to define the radial dependence of Laguerre- 
Gaussian modes one would need to perform complex am- 
plitude modulation [13], which results in the loss of op- 
tical power. 

The relative easy and the accuracy with which one 
can manipulate the azimuthal and radial dependences of 
BG modes make them a favorable choice for high-fidelity 
quantum information systems. The question is whether 
one can prepare suitable entangled states in terms of the 
BG basis. We already know that SPDC generates states 
that are entangled in terms of their azimuthal indices, 
but how does the choice of the radial dependence of the 
modal basis affect the coefficients in the OAM expan- 
sion of these states? It is preferable that these coeffi- 
cients have the same magnitude, because that represents 
a quantum state that is maximally entangled. As a re- 
sult one prefers a broad flat spectrum in terms of the 
azimuthal index (i.e. a large spiral bandwidth [Til]). 

In this paper we investigate the spectrum of BG modes 
generated by the SPDC process. Previously, the spi- 
ral bandwidth was calculated for the Laguerre-Gaussian 
modes 61 ] , but these calculations mostly focused only 
on the shape of the spectrum, discarding the magnitude. 
Ling et al. [lH attempted an analysis that gives the mag- 
nitude of the process, however their analysis made sev- 
eral assumptions and simplification that eventually lead 
to expressions that are of limited use. We develop the for- 
malism for this calculation from first principles (Sec. HI]) 
and then simplify the expression (Sec. lIIl] ) by considering 
certain special cases, making all assumptions explicit. In 
the end we obtain expressions that not only provide the 
shape of the spectrum, but also its magnitude. Using 
a generating function for the angular spectra of the BG 
modes (Appen. [B]). we then apply the expressions in the 
case where the modal basis is the BG modes (Sec. IIV[) 
and show that one can obtain a suitably flat spectrum for 
particular choices of the parameters. We end with some 
conclusions in Sec. |V] 
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II. TRANSITION PROBABILITY 

The OAM spectrum after SPDC is determined by the 
probability to measure a pair of down-converted pho- 
tons in a particular final state (pair of transverse spa- 
tial modes). If we assume that the quantum state after 
the crystal is given by a density operator p, then the 
probability for this density operator to contain a partic- 
ular final state is given by the trace tr{ppf} over the 
product of this density matrix with some final pure state 
PS = l^/X^/l that we are looking for. Assuming that 
the density operator evolves according to 

ihd t p(t) = [p(t),H I (t)], (i) 

in the interaction picture, where Hj(t) is the interac- 
tion Hamiltonian, and that the initial states is also pure 
we find that the trace can be ex- 



P(.h) = |*< 
pressed by 



tr {pp f } 



(tf/lfl/fa)!**) dti 



(2) 



The quantity to be considered is, therefore, the scattering 
amplitude A4, given by 



M = 



1 



to 



(Ml dtl. 



(3) 



If we view the SPDC process as a Feynman-diagram, then 
Eq. ([3]) represents the vertex. Implicit in this interaction 
is the requirement to conserve momentum and energy. 
Although similar to Fermi's golden rule, which was used 
in Ref. [16[ to compute the OAM spectrum, the expres- 
sion in Eq. (|3]) also contains an integral over time and, 
when substituted into Eq. ([2]), yields a true probability 
and not a transition rate. 



A. Decoupling the overlap integral 

The initial and final states are 'on shell' in the sense 
that the momentum and energy of these states obey the 
vacuum dispersion relation (shell condition) uj = c|k|. 
As a result one can insert identity operators, defined in 
terms of the momentum basis, 



1 = 



E / IM<M 



d 3 fc 



to (2tt) s 



(4) 



where the momentum basis elements obey the orthogo- 
nality condition [ijj 

(ki,r|k 2 ,s) =w (27r) 3 ^ s 5(k!-k 2 ), (5) 

with £(•) being the Dirac delta function. The cj-factor in 
Eq. ([5]) comes from the requirement for Lorentz invari- 
ance 

Identity operators, as defined in Eq. (QJ, are now in- 
serted for the pump, signal and idler beams, respectively. 
The initial state then becomes 



m=Y, J |k 3 ,r)* p (k 3 ,r) 



w 3 (2tt) 3 : 



(6) 



where 

* p (k 3 ,r) = (k 3 ,r|^ m ), 
and the final state becomes 

<*/l = £ / *:(ki,p)*J(k 3 ,ff) 

P,q J 



(7) 



/l 1 I d3fc l d3fc 2 , QS 

x ki,p,k 2 ,<? -— ¥ -—3, 8) 



where 

*:(ki >P )*j(k 2 , ff ) = (* / |k 1 ,p,k 2 ,?). (9) 

Substituting Eqs. ([6]) and JSJ into Eq. ([3]), we obtain 
M = E /*:(ki,p)**(k 2 ,g) V* p (k 3 ,r) 



d 3 fcx 



d 3 fe 



"wi (2tt)3 uj 2 (2tt)3 uj 3 (2tt)3' 



(10) 



where the phase space integrations share the same inte- 
gral sign and 



V 



1 r* 

- / k 1> p,k 3>9 JJ>(ti) k 3 ,r) dh, (11) 
h Jt Q 



which is interpreted as the Feynman rule for the vertex. 

In this way we have now decoupled the calculation as- 
sociated with the vertex itself from the overlap of the 
particular spatial modes of the pump, signal and idler. 
One can now evaluate the integrals for the vertex rule 
once and for all and afterwards use the same result for 
any particular set of initial and final states given in terms 
of the transverse spatial modes of the pump, signal and 
idler beams. 



B. The interaction Hamiltonian 

Assuming that one can neglect the frequency depen- 
dence of the second order nonlinear susceptibility x^ 2 \ 
the interaction Hamiltonian is given by 



Hi(t) 



Jv 



(12) 



where the integration runs over the volume V of the non- 
linear medium. The electric field vectors are given in 
terms of their components, which means that the indices 
represent the three spatial coordinates a,b,c = {x,y,z} 
and repeated indices are summed over. The expression 
for the quantized electric field is given by 
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E(x,t) = EH(x,()+EW(x,i) 

Ef n s (k)haj f . . . d 3 fc 

/ V — 2^ — Vs ^ ' a H k ) ex P[ iwt - m ^ k )( k ' x )] 7^X3 

a|(k) exp[— iuit + in s (k.)(k ■ x)] 



E 



?i s (k)fia; 1 
2e " 



d 3 k 



(27T) 



(13) 



where ?7 s (/c) represents the polarization vector for a given 
spin state, denoted by s, and a given propagation direc- 
tion, denoted by k. We made the spin state dependent re- 
fractive index n s (k) of the medium explicit by expressing 
the electric field in terms of the vacuum momentum ba- 
sis, which obeys the vacuum dispersion relation oj = c|k|. 
The creation and annihilation operators obey the follow- 
ing commutation relations 

[a.fkO.atpca)] = (27r) 3 ^ r SQu - k 2 ), (14) 

which implies that, for consistency with Eq. ([5]), 

(M = ^ (o|<z s (k) 



t(k)|0). 



(15) 



C. Simplifying the vertex rule 



When the electric field operates on one of the basis 
states defined in Eq. (|T5|) . one obtains 



E(x,i)|k',r) = \ J2 



2e 



?y s (fe)a s (k) exp(iwi — m s (k)k ■ x) + h.c. 



d 3 fc 
(2tt) 3 



w'at(k')|0) 



2eo 



r] r (k')exp(iuj't - m s (k')k' • x) w' |0) + 



(16) 



where the additional terms that we neglect are two pho- 
ton states that are orthogonal to the final state. Substi- 
tuting Eqs. (0, (H3J) and (JT5J into Eq. (JIT} and applying 
Eq. (jTHJ) , one obtains 

x J* np(ki)n 9 (k 2 )r7, r (k 3 ) wiw 2 w 3 

x exp(iAw^i — «Ak • x) d 3 x dti (17) 

where the 6 comes from all possible ways in which the 
electric field operators in the interaction Hamiltonian can 
be contracted with the initial and final states, together 
with the fact that the second order nonlinear suscepti- 
bility is symmetric with respect to any permutation of 
the indices; ki, k 2 and k 3 denote the propagation vec- 
tors associated with the signal, idler and pump bases, 
respectively; p, q and r denote the spin states of the sig- 
nal, idler and pump beams, respectively; a, b and c are 



spatial indices that are summed over as in Eq. (|12[k 

AbJ = IO3 — U)\ — UJ2 (18) 

and 

Ak = n r (k 3 )k 3 - n p (ki)ki - n 9 (k 2 )k 2 . (19) 

Only the exponential function in Eq. (jTTJ) contains de- 
pendencies on the spatial coordinates and time. The vol- 
ume V is assumed to be infinite along the transverse di- 
mensions x and y, and finite along the general propaga- 
tion direction z. Therefore 

/ exp(-iAk • x) d 3 a; = (2vr) 2 i 5(Ak x )S(Ak y ) 
Jv 

xS(Ak z L), (20) 

where 



S{Ak z L) = 




5 



and L is the thickness of the nonlinear medium along the 
z direction. 

If there are no restrictions on the time for the inter- 
action, other than that which is imposed by the spatial 
extent of the medium, we can extend the time interval to 
infinity, which gives 



/oo 
exp(iAurti) dtx = 2ttS(Aw). 
-OO 



(22) 



(Other restrictions could include the short duration of 
the laser pulse. Here we assume a CW source.) 
The expression for the vertex rule then becomes 

V = C 3 / 2 i W (k 1 ,k 2 ,k 3 )™(2 7 r) 3 

UI3 

x5(Au})5(Ak x )S(Ak y )S(Ak z L), (23) 

where g pqr is a dimensionless effective nonlinear coeffi- 
cient for the vertex, given by 



5p9I .(k 1; k 2 ,k 3 ) = x ( 2vZ*(h)ri h q *(k2)v c r(h) 1 '' 



2c 3 e 

X 3w§ y n p (k x )n q (k 2 )n r (k 3 ) (24) 

and which depends on the polarization states of the sig- 
nal, idler and pump, as denoted by the subscripts p, q 
and r, respectively. Henceforth, we'll distinguish between 
the different combinations of spin states that give phase 
matching, leading to either type I or type II phase match- 
ing. The expressions for the nonlinear coefficient g and 
Ak in these two cases are provided in Appen. [A] for the 
benefit of the reader. 

The vertex rule, as expressed in Eq. (|23l) . contains no 
more integrals. It does contain a number of Dirac delta 
functions that impose energy and momentum conserva- 
tion (to some extent) and a sinc-function S(Ak z L) that 
imposes the phase matching condition. 



D. General expression 

We now substitute V, given in Eq. (|23|) . into A4, given 
in Eq. (fTU|) , to obtain 

M = c^ 2 L J 5 (k 1 ,k 2l k 3 )^(k 1 )^(k 2 )^(k 3 ) 
x (2Tif 5 -^^-5{Ak x )8{Ak y )S{Ak z L) 



d 3 fci d 3 fc 2 d 3 fe 3 
J2nf J2^f (^O 3 ' 



(25) 



where we dropped the explicit spin state dependencies 
and the summation over all the spin states, because we 
either consider type I or type II separately, so that g can 
either be g ooe or g eoe as given in Appen. [A] For type I 
phase matching 

Ak z = n e ff(6 x ,U3,03,<p3)k Z 3 - n (u\)k z \ - n (w 2 )/c z2 , 

(26) 



and for type II phase matching 

Ak z = n ef f(d x ,uJ3,0 3 ,(f>3)k z3 - n e ff(9 x ^i,0 1 ,(j) 1 )k zl 



-n (w 2 )fc z2 . 



(27) 



The expression in Eq. (|25l) is quite general, but it may 
not be convenient for calculations. Therefore, we subse- 
quently introduce a number of simplifications. 



III. SIMPLIFICATIONS 

A. Choosing a propagation direction 

Due to the vacuum dispersion relation, any three of 
the four quantities lo, k x , k y and k z will fix the remain- 
ing quantity. Therefore, the phase space integrals for 
propagating (on shell) fields only run over three of these 
quantities, usually k x , k y and k z . For an optical beam, 
continuously propagating through a linear medium, it is 
more convenient to specify lo, k x and k y and thereby fix 
k z . The reason is that the system's evolution is consid- 
ered as a function of the propagation distance and not as 
a function of time. The input is therefore specified for all 
time as a two-dimensional function on a plane perpendic- 
ular to the propagation direction. Such an input can also 
be expressed in the Fourier domain as a two-dimensional 
function of k x and k y for a particular lo. (This implies 
a monochromatic approximation, however one can also 
specify the input as a function of lo.) In other words, we 
prefer to integrate over lo instead of k z . The orthogonal 
basis for the quantum states is then defined by |K,w, s) 
instead of Ik, s). where 



K = k x x + k y y 



(28) 
(29) 



and s denotes the spin state. By fixing a specific direc- 
tion for propagation (the z-direction) we explicitly break 
rotation invariance and, by implication, also Lorentz in- 
variance. 

From the vacuum dispersion relation it follows that 



dfc. 



c 2 k. 



dco. 



(30) 



In the applications that we consider the angular spectrum 
of the beam only contain nonzero components in the pos- 
itive k z region. Since both sides of the to axis map into 
the positive side of the fc z -axis the change of integration 
variables in Eq. (|30|) is well defined. Applying this change 
of variables to Eq. (|25|). we obtain 



M 



C 3 /2 L J $*( Kl) Wl )*J(K 2l w 2 )* p (K 3 , wa) 



X5(ki,k 2 ,k 3 



UJlU) 2 

i 

UJ3 



x (2Tif5{Aio)5(Ak x )5{Ak y )S(Ak z L) 
d 2 fci Au>i d 2 k 2 duj 2 d 2 k 3 dw 3 



c 2 fc zl (2^) 3 C 2^ 2 (2^) 3 C 2 fc 23 (27T) 3 



(31) 
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where 



i z2 



';3 



UJ 2 

c 2 


~ ^xl 


- k 2 
K yl 


J 2 

c 2 


~ ^x2 


- k 2 
K y2 


UJ 2 

c 2 


~ ^x3 


- k 2 
K y3 



-n (w 2 )V -f - k 2 x2 - k 2 



vi- 



and for type II phase matching 
Ak z = n eff (0 x ,uJ3,O 3 ,(f> 3 ) 



(32) 



For type I phase matching 



Ak z = n eff (6 x ,u; 3 ,9 3 ,cf> 3 )\l - k* 3 - k% 3 



(33) 



k 2 - k 2 



„ 2 ^x3 ^yZ 



-n ef f(6x,ui,0i,<t>i) ~ k 2 xl - k 2 



-n {u2)\l^-k 2 x2 -k 2 y2 . 



(34) 



The normalization now reduces to 

|2 d 2 fc do; f 9 d 2 fc 



'^ (K '" )r c 2 fc z ( 27 r)3 
where we use the fact that 



ic(K)r 



(27T) 



\H(lj - uj p ;Su})\ 2 = 27r. 

Z7T 



(38) 
(39) 



Substituting Eq. (|3"5j) into Eq. (|3"TT) and evaluating the 
integral over ui 3 , we obtain 

X = c 3 / 2 L / *:(K 1 ,^ 1 )*.*(K 2 , W2 )G(K 3 ) 



xg(ki, k 2 , k 3 )H(oj! +uj 2 - w p ; <M 

x<5(A/c a; )(5(Afcj / )S'(Afc ZJ L) 

d 2 fc 3 d 2 fci dwi d 2 fc 2 dw 2 



WiW 2 

Wl + w 2 



cVk^ c 2 fc zl (27r)3 c 2 fe 2 (2^)3 



(40) 



where Ak z is given by Eqs. (|3"3")l and (jM)) with the re- 
placement W3 — > uj\ + cj 2 , and 



%3 



(Wl + W 2 ) 5 



£- 2 — k 2 
ft a:3 S3' 



(41) 



B. Monochromatic pump 

We now assume that the pump is monochromatic. For 
this purpose we define the momentum space wave func- 
tion of the pump by 



Vv( K > w ) 



k z G{¥L)H(w-u p ;5w). 



(35) 



where uj v and Su> are, respectively, the center frequency 
and the (small) bandwidth of the pump laser. (One can- 
not define H(ui) as a Dirac delta function, because that 
would lead to a squared Dirac delta function under the 
integral in the normalization condition, which gives a di- 
vergent result.) To satisfy the normalization requirement 
for the momentum space wave function 



1 / nt m 2 d2fc duj 1 
'^ (K ' w)l ?M2^ ' 



we define 



H(u;Su) 



2VV/4 



■ exp 



25a; 2 



(36) 



(37) 



The frequency spectrum in Eq. Q37p does not actually 
enforce the pump to be monochromatic, unless we take 
6u -> 0. 



C. Degenerate, collinear type I phase matching 

At this point we restrict ourselves to type I phase 
matching, with degenerate signal and idler frequencies 
(ui = L02 = = uip/2) and collinear signal and idler 
beams. Due to the phase matching condition for collinear 
beams, the effective refractive index for the pump beam 
must be equal to the ordinary index at the degenerate 
frequency [n e f / (Ox, w p , 0,0) = n (uJd)}- We'll therefore 
denote all refractive indices simply by n . Moreover, as- 
suming that the nonlinear coefficient is a slow varying 
function of the three propagation vectors and that the 
beams are paraxial, we replace g(ki, k 2 , k 3 ) with a con- 
stant g and pull it out of the integral. 

The degeneracy condition is implemented by replacing 



*J(K a ,W2) 



ft,iM*(Ki)ff(ui -u> d ]Suif) 
~k^M*(K 2 )H(u> 2 -u d ;5LJ f ) (42) 



in Eq. (|40|) . where H(-) is defined in Eq. (|37|) and Sujj is 
the bandwidth of the line filter that imposes degeneracy. 
Assuming that the bandwidths are small enough one can 
simply substitute tjj\ — cj 2 = ujd everywhere, except in 
the iJ-functions. Evaluating the integrals over uj\ and 
w 2 , one obtains 
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H(u)i - uj d ; 6uf)H(u 2 - 6b}f)H(u3i + uj 2 - 2uj d ; Slu) ^p- = 

Z7T Z7T 7T ' 



/5w 



<5o; 2 



2(5w 2 



(43) 



where the approximation follows from the fact that the 
line filter bandwidth is usually much larger than the 
pump bandwidth 8u3 -C <$wj. 

In the paraxial limit the remaining /c z -factors become 



(k z ik z2 k z3 



,-1/2 



3/2 



[l + O(0 2 )] 



(44) 



where 6(= \ p /irdo) is the pump beam angle, with do 
being the radius of the pump beam profile at its waist. 
Therefore, we use only the leading order term for this 
factor. 

The degenerate, collinear type I phase matching con- 
dition implies that 



6{Ak x )6(Ak y ) 



1 



S(K 3 — Ki K 2 



(45) 



Using Eqs. ([42|). ([43]). (1441) and (|45|l. and evaluating the 
phase space integrals for the pump to eliminate the re- 
maining Dirac deltas in Eq. ([40]) . we obtain 



M 



Lg I 5lj 




:G(Ki +K 2 )5(Afc z L) 



m;(K!)m;(k 2 ) 

d 2 fc x d 2 fc 2 



(2tt) 2 (2tt) 2 



(46) 



where Ak 7 



n a (k 



z3 



k z i - k zl ) and 



'• 1 — »' c 2 ""art Si 



h 7 



d _ u2 _ L.2 
c 2 K x2 K y2 



Au 2 

/,-., = \/^-(fc*i + M 2 -(^i + M 2 - (47) 
In the paraxial limit Ak z simplifies to 

(hci - k x2 ) 2 + (k y i - k y2 ) 2 



n D c- 



4w rf 



(48) 



It is often more convenient to do the calculation in terms 
of spatial frequencies instead of the propagation vector 
components. Therefore, we define the propagation vector 
in terms of spatial frequencies by 



K = 2ir(ax + by). 



Then Ak z becomes 
Ak z = — 



[(oi - a 2 ) 2 + (bi - h) 2 ] , 



(49) 



(50) 



where 



Hence 



K = TTTloXpL. 



S(Ak z L) = 



sin {2k [( ai - q 2 ) 2 + (h - b 2 ) 2 ] } 
2K[(ai-a 2 ) 2 + (6i-6 2 ) 2 ] ' 



(51) 



(52) 



IV. BESSEL-GAUSS OAM SPECTRUM 

We now evaluate the integrals in Eq. (|4"6"|) for a par- 
ticular choice of momentum space wave functions for the 
pump, signal and idler beams, respectively. It is assumed 
that the pump beam has a Gaussian profile, so that 



G(K) = V2^d exp [■ 



■b 2 )] 



(53) 



where do is the beam waist radius. The momentum space 
wave function in Eq. (|53|) satisfies the normalization con- 
dition given in Eq. ([38]). 

The signal and idler beams are assumed to be BG 
modes. For convenience we'll use the generating function 
for the angular spectra of the BG modes (see Appen. [5J , 
which is given by 



G{a,b) 



2ird n exp{-ir 2 d 2 n (a 2 + b 2 ) - d 2 n h 2 n /A 
irhnd^ [cos(a„)6 - sin(a„)a] } , (54) 



where a and b are the spatial frequency coordinates, h n 
is the radial mode parameter of the BG modes, d n is the 
radius of the Gaussian envelope and a n is the generating 
parameter which is used, according to Eq. (|B3I) . to gener- 
ate the angular spectrum of a BG mode with a particular 
azimuthal index £. The subscript n(= 1, 2) represents the 
signal or idler beam. We set z = inside the generat- 
ing function, thereby assuming that the waists of the BG 
modes are located in the nonlinear crystal where these 
modes are generated. 
Substituting 

m;(K!)m;(K 2 ) = e s *(ai,&i)£*(a 2 ,& 2 ), (55) 

and Eqs. flS^) and (53) into Eq. (|35|). and using Eq. ([4"9"]l. 
we obtain the integral expression for the generating func- 
tion of the scattering amplitudes 



A = 



Qo J exp {-d 2 n 2 [(oi + a 2 ) 2 + (6 a + b 2 ) 2 } } 
sin {2k [(ai-a 2 ) 2 + (fei-6 2 ) 2 ]} 



2k [(ai 



ai) 



(bi 



xG*(ai,b 1 )Gt(a 2 ,b 2 ) dai d&i da 2 d& 2 , (56) 
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where 



(57) 



A. Result after integration 



be done for both the signal and idler beams, the relevant 
expression is 



Mr 



1 



4tt 2 



exp(mai) exp(ima2) 

io Jo 
xA(ai — Oi2) dai da2. 



(64) 



To expedite the evaluation of the integrals in Eq. (1551) , 
we'll assume that h\ = h% = h, and also that d\ = d,2 = 
loq. Moreover, we assume that the pump wavelength X p 
is much smaller than any of the other dimension param- 
eters. Under these conditions the integrals evaluate to 



A 



n 0l 2 



■ exp 



l + 2 7 2 



irdfa(l + 2 7 2 )£ 2 cr 2 
x - (1 + irj)t;a] exp[(l + ir])£a] 
- (1 - irj)&] exp[(l - ir))£a]} 



where 



7 

n 



do_ 

h 2 oj 2 
n X p L 

1 . (ax- 
— sin 

2 V 2 



0:2 



and 



Lg 



(58) 

(59) 
(60) 
(61) 

(62) 
(63) 



The overall absolute magnitude of the probability for 
producing down-converted photon pairs is represented by 
the dimensionless quantity Q\ given in Eq. (|63[) . Apart 
from some numerical constants and the nonlinear coef- 
ficient g, it depends on the ratio of the length of the 
nonlinear crystal L and the radius of the pump beam do, 
as well as on the square-root of the fraction of the band- 
width of the pump compared to the center frequency of 
the pump. 



B. Extraction of the different orders 

The coefficients for the different orders are extracted 
with the aid of Eq. (|B3|) . Since the extraction needs to 



However, since A only depends of the difference in angle 
ol\ — a,2 and since these angles are cyclic variables, one 
can redefine one of the angles as U2 — ol\ — 03, so that 
A becomes independent of a\. The integral over a± then 
gives a Kronecker delta function, 



1 f 27T 

— / exp[— i(n + m)ai] d«i = <L 
2tt Jo 



(65) 



which shows that orbital angular momentum is conserved 
in the SPDC process. The remaining integral over 0:3 is 
given by 



M m = — A(a 3 ) eyLp(-ima 3 ) da 3 . (66) 
2?r Jo 




FIG. 3: The curves for the normalized magnitudes of the coef- 
ficients l-MmpOl as a function of X for a number of different 
coefficients (m — 0...5). 



Upon evaluating the integral over 0:3, we obtain 



Mm = (-irfw P (-f){ T ^/ m ^ 



24r 



[(1 + 2 7 2 )m 2 + (1 + 2 7 2 )(m + X) 2 + 4X~f 2 + 2(2 7 2 - l)m]J„ 



(67) 



up to second order in rj, where I m (-) is the modified Bcssel 
function and 



X = 



£7 2 



"o u o 



l + 2 7 2 



2d 2 ' 



(68) 



Ignoring the ?7 2 -terms in Eq. (1671) . we see that the domi- 
nent behavior of the coefficients is given by 



\m,,a vi| • i-yl^'fy 



(69) 




FIG. 4: The curves for the OAM bandwidth as defined in 
Eq. (f70|) as a function of X. 

In Fig.[3]we show the curves for Eq. ([69)) as a function 
of X for a number of different coefficients. For X = 



the first coefficient (m = 0) is 1 and the rest are zero. 
As X becomes larger the magnitudes of the coefficients 
become closer to each other. It is therefore desirable to 
make X as large as possible. 



C. OAM Bandwidth 

One can compute the OAM bandwidth of the down- 
converted beams by exploiting the generating function in 
Eq. <J5SJ). We define the OAM bandwidth by 



B 



E m 2 \M r 
E\M m \ 



1/2 



(70) 



Formally the generating function for the coefficients can 
be written as 



(71) 



The conservation of angular momentum of Eq. (|65[) turns 
the double summation into a single summation. 

The bandwidth as defined in Eq. (f70| can be computed 
directly from the generating function in Eq. (|58[) as fol- 
lows 



B = 



J\d a A(a)\ 2 da 
, f\A(a)\ 2 da , 



1/2 



(72) 



The resulting expression for the bandwidth, up to sec- 
ond order in 77, is 



B = - 



1 


~xh{xy 


1/2 | 


2 


[ MX) _ 





-(4yx + 2X + ioy* + 1 



?7 2 (l + 2 7 2 ) 
96 7 4 
XI (X) 



2(1 + 2 7 2 ) + (4 7 2 A + 2AT + 6 7 2 - 1) 



XhjX) 
Io(X) 



h{X) 



(73) 



where X is given in Eq. f|68[) . The leading order term 
gives the bandwidth purely as a function of X. In Fig. 0] 
we show the curves for the OAM bandwidth as defined in 
Eq. (l7Cfl) as a function of X. We see that the bandwidth is 
a monotonically increasing function of X, which confirms 
the observation at the end of the previous section. Since 
X is proportional to the square of the transverse scale 
parameter for the Bessel functions h, we find that a large 
beandwidth is obtained for large values of h. 



V. CONCLUSIONS 

The OAM spectrum for BG modes generated by the 
SPDC process is calculated. This calculation is preceeded 
by a general derivation of expressions for the spectrum 
generated by the SPDC process in any transverse spa- 
tial modal basis. Through a series of simplifications, we 
eventually arrive at expressions for the OAM spectrum 
for BG modes in degenerate, collinear SPDC with type 
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I phase matching. The calculation of the OAM spec- 
trum is done using generating functions for the angular 
spectra of the BG modes. We obtained a closed form ex- 
pression for the bandwidth of this OAM spectrum. The 
results show that the magnitudes of the coefficients in 
the BG modal expansion of the quantum states of the 
down-converted photon in terms of become close to be- 
ing equal, provided that the transverse scale parameter 
for the Bessel functions is as large as possible. 

The expressions for the OAM spectrum also provide 
the overall absolute magnitude of the spectrum, which is 
represented by the dimensionless quantity 51 1 in Eq. ([55]). 
however we did not investigate this quantity any further. 

The next step would be to perform physical experi- 
ments in which the coefficients for the different sets of 
modes are measured. In this way one can check whether 
the shape and the absolute magnitude of the observed 
spectra agree with the theoretical results of this paper. 
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Appendix A: Effective nonlinear coefficients 

Considering the case of a negative uniaxial crystal for 
which extraordinary index n e is smaller than the ordinary 
index n , one obtains the following two phase matching 
conditions: 



• Type I phase matching, for which 

9ooe = \J n e ff(0x,U3, #3, </>3)rc (wi)n (w 2 ) 

x 3w 3 2 ^- X ^X* (ki )rfc (« {k) (Al) 

and 

Ak/ = n e// (6>x,w 3 , 03,03^3 

-n (wi)ki - n (cj 2 )k 2 (A2) 

• Type II phase matching, for which 

9eoe = \J n e ff{9x,U3, 6> 3 , (j> 3 )n e ff(6 x , ui, 0~i, 4>i)n (ui 2 ) 

x ^fy^x&C Cki)v b * (h)v c e (k) (A3) 

and 

Ak// = n e// (6> x ,w 3 ,6l3,03)k3 

-n e ff(9x,ui,9i,<t>i)ki - n (w 2 )k 2 . (A4) 

In the above expressions we define the propagation vec- 
tors by 

k = k (x sin cos <j> + y sin sin <j> + z cos 0) (A5) 
and the optic axis by 

a = y sin Ox + z cos 6x ■ ( A6) 

The effective refractive index n e ff(-) is given as 



n e ff{0x,u,t 



n e (w)n (w) 



\J n 2 (uj) + [n 2 (w) — n 2 (w)] (sin Ox sin sin <fi + cos Ox cos 0) 2 



(A7) 



The parameters with the subscripts 1, 2 and 3 in 
Eqs. (|A1I) - (|A4I) represent the parameters for the signal, 
idler and pump beams, respectively. 



Appendix B: Generating function for the 
Bessel-Gaussian modes 

In normalized coordinates the Bessel-Gaussian (BG) 
modes are given by 

M? G (r,<M;x) = sj^h (yzu) ex P^ 



where the normalized radial coordinate is given by r — 
(x 2 + y 2 ) 1 / 2 /loq, cf> is the azimuthal coordinate and the 
normalized propagation distance is t = z/zr] t is the 
mode index (a signed integer) and x(= huio) is a normal- 
ized parameter for the size of the Bessel mode. The ini- 
tial radius of the Gaussian profile is ujq and the Rayleigh 
range is Zr — ituj 2 /\. There is an additional propaga- 
tion phase factor associated with every BG mode given 
by exp(— ibz), which is not shown in Eq. (|B1|) . The 
wavenumber k and the parameters b and h are related 
by k 2 — b 2 + h 2 . So b and h are, respectively, the mag- 
nitudes of the longitudinal and transverse components of 
the dominant wave vectors in the spectrum of the BG 
mode. 

The generating function for the BG modes for a par- 
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ticular \ is given by 



modes and is given by 



77 1 — 



— exp 



iAxrsin((f) — a) — 4r 2 + ix 2 t 
4(1 - it) 



(B2) 

where a is an angular generating parameter that is used 
to generate a BG mode with a specific value of the az- 
imuthal index I. To generate such a mode one performs 
the following calculation 



Mf G 



1 

2t7 



2tt 



G exp(ito) da. 



(B3) 



The Fourier transform of the generating function gives 
the generating function for the angular spectra of the BG 



/oo 
Q exp[i27r(au + bv)] da db 
-oo 



'2n exp 



-7T 2 /z 2 (l — it) + TrxfJ' sin(V — a) — 



(B4) 



where fi and v are, respectively, the normalized radial 
coordinate and the angular coordinate for the spatial fre- 
quency domain. The angular spectra of the BG modes 
are generated from Eq. (|B4[) in the same way as in 
Eq. (55) . 
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